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ABSTRACT 



Recently, peculiar velocity measurements became available for a new sample 
of galaxy clusters, hereafter the SCI sample. From an accurately calibrated 
Tully-Fisher relation for spiral galaxies, we compute the rms peculiar velocity, 
V rms , and compare it to the linear theory predictions of COBE-normalized 
low-density and open CDM models (ACDM and OCDM, respectively). 
Confidence levels for model rejection are estimated using a Monte Carlo 
procedure to generate for each model a large ensemble of artificial data sets. 
Following Zaroubi et al. (1997), we express our results in terms of constraints on 
the (Q ,n pr , h) parameter space. Such constraints turn into a^o' 6 = 0.50lgj4 
at the 90% c.L, thus in agreement with results from cluster abundance. We 
show that our constraints are also consistent with those implied by the shape 
of the galaxy power spectrum within a rather wide range for the values of the 
model parameters. Finally, we point out that our findings disagree at about 
the 3<t level with respect to those by Zaroubi et al. (1997), based on the Mark 
III catalogue, which tend to prefer larger Vt Q values within the CDM class of 
models. 
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1. Introduction 

Peculiar velocities of clusters have recently been used by several authors to set stringent 
constraints on cosmological models (e.g., Croft & Efstathiou 1994; Cen, Bahcall & Gramann 
1994; Bahcall & Oh 1996; Moscardini et al. 1996). Although clusters sample the large-scale 
flows much more sparsely than galaxies, their peculiar velocity can be measured more 
accurately if distances are available for several cluster galaxies. Several of such previous 
analyses were, however, based on non-homogeneous compilations, with cluster velocities 
taken from different parent samples. 

In this Letter we analyze the new sample of cluster peculiar velocities described 
by Giovanelli et al. (1997a; SCI hereafter), consisting of accurate and uniform I-band 
photometry and velocity width measurements for about 800 spiral galaxies in the fields of 24 
clusters. Of the 24 clusters, we consider 18, pruning the six paired clusters (A2197/A2199, 
S805=Pavo II/Pavo, and A2634/A2666) in order to avoid possible ambiguities in 
membership assignment (see Giovanelli et al. 1997a). 

Peculiar velocities for these clusters were determined from a a well-calibrated I-band 
Tully-Fisher relation constructed from all 24 clusters as described in Giovanelli et al. 
(1996b). An earlier version of this sample has been already analyzed by Bahcall & Oh 
(1996) and Moscardini et al. (1996), who compared the results with numerical simulations 
of several CDM-like models. Within this class of models, both analyses consistently 
showed that this data set favors a low-density Universe, with 0.2^ f2 ~ 0.4. Furthermore, 
Moscardini et al. (1996) also compared the SCI sample with Hudson's (1994) compilation of 
cluster velocities. They found that the SCI sample provides systematically smaller velocities 
than those of Hudson, again suggesting a low value of VLq. 

The analysis that we present in this Letter is entirely based on linear theory, the 
reliability of which to describe cluster motions is briefly discussed. Model predictions are 
worked out for purely CDM models with Q < 1 and both flat and open geometry. Avoiding 
the need to resort to numerical simulations allows us to probe the model parameter space 
in a much more accurate way. The resulting constraints on the CDM models are also 
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compared with those derived from the cluster abundances, the galaxy power-spectrum 
shape and the Mark III data as analyzed by Zaroubi et al. (1997) 



2. The analysis 



Our analysis is based on comparing the rms cluster velocity, V rms from the SCI sample 
and for models. Its observational estimate for the SCI clusters gives a one dimensional 
K-ms = 266 ±30 km s -1 , where the uncertainty represents the la scatter over 10 5 Montecarlo 
realizations of the real sample, each one generated from an a priori Gaussian distribution 
having the same V rms as the SCI data set and velocities convolved with the observational 
errors. 



Linear gravitational instability predicts the one-dimensional rms velocity to be 
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where /(f^o) — 6 ' * s the m °del power spectrum and W(kR) is the window function 
that specifies the "shape" and the size R of the linear density fluctuations, which generate 
clusters. Bahcall, Gramann & Cen (1994) and Croft & Efstathiou (1995) verified that eq.fll]) 
provides a rather good fit to the cluster rms velocity generated by N-body simulations. 
By using the Gaussian window W(kR) = exp(—k 2 R 2 /2), we found that eq.([l|) provides 
the best fit to the Borgani et al. (1997) N-body outputs for a variety of models by taking 
R = 3.9/i -1 Mpc (Hq = lOO/ikms^ 1 Mpc -1 ), corresponding to a typical cluster mass 
M c i ~ 2.6 x lO 14 r2o^- _1 ^0- We adopt this value in the following analysis. 

We take the power spectrum to be P(k) = Ak Upr T 2 (k), where 

T(q) = ln(1 2 + 3 ^ 34g) [l + 3.89g+(16.1g) 2 + (5.46g) 3 +(6.71g) 4 ]" 1/4 (2) 

is the CDM transfer function provided by Bardeen et al. (1986). Here q = k/Th where 
T = £l hexp(— fi;, — (2/i) 1 / 2 fife/fi ) is the "shape" parameter, which takes into account 
the presence of a non-negligible baryon fraction, Qf, (e.g., Sugiyama 1995). We take 
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VL b h 2 = 0.024 (e.g., Tytler et al. 1995). The constant A is fixed by the 4-year COBE 
normalization recipes of Bunn & White (1996; see their eqs. 29 and 31) and Hu & White 
(1997; see their eq. 6) for flat low-density (ACDM) and open (OCDM) models for both 
vanishing and non-vanishing tensor mode contributions to CMB anisotropies. In the 
case of ACDM, Bunn & White (1996) consider the case T/S = 7(1 — n pr ) for the ratio 
betweem the quadrupole moments of the tensor and scalar modes in the expansion of the 
angular temperature fluctuations as generated by power-law inflation (e.g., Crittenden et 
al. 1993, and references therein). As for OCDM, the normalization for the T/S ^ case is 
provided by Hu & White (1997) in the case of "minimal" tensor anisotropies. The density 
parameter Qq and the primordial spectral index n pr are varied within the ranges where the 
normalization fits are reliable, namely 0.2 < Q < 1 and 0.7 < n pr < 1.2 (0.7 < n pr < 1) for 
the cases without (with) a tensor mode contribution. 

The family of models we consider is specified by the three parameters (Qo,n pr , h). 
Results will be presented by keeping one of them fixed, in the form of slices of the 
three-dimensional parameter space. Colberg et al. (1997) have recently pointed out 
that cluster peculiar velocities are almost independent of the density parameter, once 
cosmological models are normalized to reproduce the cluster abundance, according to the 
recipe by Eke et al. (1996). Here we follow the different approach of imposing the COBE 
normalization, since we regard CMB temperature anisotropies as a more stable and robust 
constraint than cluster abundance for a fixed choice of model parameters. 

In order to establish the confidence level for the validity of a given model, we adopt 
the following procedure. Let Vi and <7, be the velocity and its error for the i-th real cluster 
(i = 1, . . . , 18). For a given model, we generate Montecarlo samples, each containing 
18 velocities, Vi, drawn from a Gaussian distribution, having dispersion given by eq. ([]]). 
For each sample we convolve the z-th model velocity with the observational error Oi and 
estimate the resulting rms velocity. For each sample, every cluster's velocity is newly 
estimated as a Gaussian deviate of the mean Vi and dispersion <jj, and the rms velocity of 
the sample then computed. 

For each model we generate N = 10 4 samples and then compute the fraction T of them 
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with Vj ms (j = 1, . . . , N) at least as discrepant as V°^ s with respect to their average value, 
N~ l Yjj Vrms'- the smaller the value of J 7 , the smaller the probability that V°^ s is generated 
by chance by that model, the larger the probability V = 1 — T that the model itself is 
rejected. 

Figure 1 shows the results of our analysis for scale-free (i.e., n pr = 1) ACDM and 
OCDM models, also comparing them to other observational constraints. The contours 
indicate the iso-probability levels for the model exclusion. The outermost contour is for 
V = 90% confidence level, while different levels are equi-spaced in logarithmic units by 
A(\ogV) =0.1. The heavily shaded area indicates the la confidence level from the cluster 
abundance by Eke et al. (1996). By fitting the X-ray cluster temperature function with 
CDM model predictions, they found <r 8 fig = 0.52 ± 0.04, with a = 0.52 - 0.13fi for ACDM 
and a = 0.46 — 0.10f2o f° r OCDM (erg is the linear rms density fluctuation within a top-hat 
sphere of 8/i -1 Mpc). Q The medium-weight shaded area is the 95% confidence level 
from the fitting by Liddle et al. (1996) to the shape of the APM galaxy power spectrum 
by Peacock & Dodds (1994): T = 0.23 - 0.28(1 - l/n pr ) with errors of about 16%. A 
consistent result has also been found by Borgani et al. (1997) from the analysis of the 
cluster distribution. The lightly shaded area shows the 90% confidence level by Zaroubi 
et al. (1997; Z97 hereafter) from the likelihood analysis of the Mark III galaxy peculiar 
velocities (Willick et al. 1996), whose results are reported in Table 1. The dashed curves 
are for different values for the age of the Universe: to = 9, 11, 13, 15, 17 Gyrs from upper to 
lower curves. 

Our results differ with respect to those by Z97. The difference is larger for ACDM 
models, for which the discrepancy is at ~ 3a level (note that the corresponding 90% 
confidence level are at most marginally overlapping), the latter favoring larger Qq values 

2 Pen (1996) recently pointed out that the scaling by Eke et al. (1996) somewhat 
underestimates the value of as at small Qq values. In particular, he claimed that for a 
ACDM model with ~ 0.35 a$ should be ~ 17% larger. We checked that the central value 
for the h interval corresponding to such an Qq increase from h = 0.66 to h = 0.74. 
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(at a fixed h). This result points in the same direction as that found by Moscardini et al. 
(1996). On the other hand, the constraints we set on CDM models are quite consistent with 
those coming from the P(k) shape and the cluster abundance, in a rather broad range of Q 
and h values. For instance, if we impose ages in the range 13^ t ~ 15 Gyrs, ACDM models 
require 0.35^ Vl < 0.50 with 0.50^ h< 0.65, while OCDM models require 0.50^ n < 0.70 
with 0.45^ h^, 0.60. On the contrary, the results by Z97 are rather discrepant with both 
such constraints, especially with the cluster abundance. 

A similar picture also emerges as tilted (i.e., n pr ^ 1) models are considered. Figure 
2 is analogous to Fig.l , but with results plotted in the n pr -Q plane, taking h = 0.65 
and 0.55 for ACDM and OCDM models, respectively. For both classes of models, taking 
T/S has the effect of (a) narrowing the permitted region in the parameter space and 
(b) decreasing the need for a tilt (cf. also Z97). Tilting P(k) breaks the degeneracy of 
the P{k) shape with the other constraints. Fixing h = 0.65 (e.g. Giovanelli et al. 1997c) 
and t — 13 Gyrs for ACDM would require f2 = 0.43; this turns into 0.85^ n pr ^ 0.95 and 
0.90^ n pr ^> 0.96 for T/S = and T/S = 7(1 — n pr ), respectively. Consistency between 
the P(k) shape and Z97 are attained for n pr ^ 1 and fi ~ 0.5, while the cluster abundance 
is still largely missed. As for OCDM models, taking h = 0.55 and t ~ 13 Gyrs implies 
Q ~ 0.65 and 0.84< n pr < 0.94 (0.84< n pr < 0.94) for T/S = (^ 0). A substantially 
larger h value would turn into too small Q values, unless t < 13 Gyrs. Again, the SCI 
cluster velocities are consistent in all the cases with the other two constraints for reasonable 
values of the model parameters. 

In order to better quantify the difference with respect to the constraints provided by 
the Z97 analysis, we fit the same combination of parameters, fi ^5o n p r = C (^50 — 2 h: 
Hubble constant in units of 50kms _1 Mpc _1 ), considered in that paper and the results are 
in Table 1. Although the shape of the relation (i.e., the values of /x and v) is quite similar, 
its amplitude C is significantly different. This confirms that, for fixed h and n pr values, 
our results favor a lower density parameter. We have also computed the best fit to the 
quantity agQ^ 6 , which fixes the amplitude of the velocity field in linear theory. We find 
og^o 6 = 0.50to]if (errors correspond to 90% confidence level), which again agrees with the 
constraints from the cluster abundance, and is significantly smaller than <t 8 ^o' 6 = 0.88 ±0.15 
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that derived by Z97. Willick et al. (1996) recently compared the Mark III data with 
velocity and density fields reconstructed from the 1.2 Jy IRAS survey. Quite interesetingly, 
they obtained a 8 ^Q 6 = 0.34 ± 0.05 (cf. their Fig. 20), thus at variance with respect to the 
results by Z96, also based on the Mark III sample, and rather consistent with our results. 

3. Conclusions 

We have performed a detailed comparison of the cluster peculiar velocities in the 
SCI catalog with those predicted by COBE-normalized CDM models, using linear theory. 
This comparison has been made by computing the rms cluster velocity, V rms , for data and 
models, and estimating the likelihood that the observed value V rms = 266 ± 30kms _1 is 
consistent with a given model. 

Confidence levels for rejecting models were determined using a Monte Carlo procedure 
which generates a large number (10 4 ) of mock samples from each model. The main goal 
of our analysis has been to impose constraints on the space of (Q ,n pr , h) parameters for 
CDM models. We have compared our results with those of Z97, and with the constraints 
that have been established from the properties of clustering of galaxies as expressed by the 
shape of the power-spectrum and recent determinations of cluster abundance. 

Our results can be summarized as follows: 

(a) Velocities of SCI clusters point toward a low-normalization model, characterized by 
cr 8 f2Q 6 = 0.50j;oj7. This result agrees with the independent constraint coming from the 
abundance of galaxy clusters. 

(b) Our results disagree at about the 3a level with those of Z97, based on Mark III, the 
latter generally indicating higher velocities and, therefore, favoring larger Qq values for 
fixed h and n pr parameters (cf. Table 1). On the other hand, we are quite consistent 
with the analysis by Willick et al. (1996), which is also based on the Mark III sample. 

(c) The results agree well with those from the analysis of field spirals in the new SFI 
sample (da Costa et al. 1997; Freudling et al. 1997). 
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The conclusions that we draw in this Letter about the values of the model parameters 
strictly hold, only for the CDM class of models. For instance, Cold+Hot DM models, 
are characterized by different power spectrum shapes and smaller COBE-normalized <j 8 
values for a fixed choice of (f2 ,n pr , h), depending on the amount and the nature of the hot 
component (e.g., Primack 1996, and references therein). We postpone to a forthcoming 
paper the analysis of a wider class of cosmological models, as well as the comparison with 
other data sets for galaxy and cluster peculiar velocities. 
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Figure captions 



Figure 1. Constraints from the SC cluster sample on the (h,Qo) plane for scale-free 
(n pr = 1) models. The countours are the levels at equal probability V for model 
rejection. The most external level corresponds to V — 90% and the spacing corresponds to 
A(log"P) = 0.2. The heavely shaded area is the constraint from cluster abundance (after 
Eke et al. 1996), the medium-weight shaded area is for the shape of the APM galaxy 
power spectrum (after Viana et al. 1996) and the lightly shaded area is from the analysis of 
Mark III velocities by Zaroubi et al. (1996). Dashed curves indicate different ages for the 
Universe: to = 9, 11, 13, 15, 17 from upper to lower curves. 

Figure 2. It is analogous to Fig. 1, but on the (n pr ,f2 ) plane. Different ages of the 
Universe are now indicated with the vertical dashed lines. 
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Table 1: Values of thw fitting parameters for the relation VL^h^n^. = C (h 50 — 2h: Hubble 
constant in units of 50kms _1 Mpc _1 ), from our analysis and from that by Zaroubi et al. 
(1997; Z97). 
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